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Abstract 

We discuss the relation of the specific heat, the energy density and the 
thermodynamic Casimir effect in the case of thin films in the three dimen- 
sional XY universality class. The finite size scaling function 9(x) of the 
thermodynamic Casimir force can be expressed in terms of the scaling func- 
tions h'{x) and h{x) of the excess energy density and the excess free energy 
*73 ■ density. A priori these quantities depend on the reduced temperature t and 

the thickness Lq of the film. However finite size scaling theory predicts that 
the scaling functions depend only on the combination x = tfLo/^o] 1 ^; where 
v is the critical exponent and £o the amplitude of the correlation length. We 
exploit this fact to compute 9 from Monte Carlo data for the excess energy 
' density of the improved two-component (f) A model on the simple cubic lattice 

with free boundary conditions in the short direction. We repeat this exercise 
00 , using experimental data for the excess specific heat of 4 He films. The finite 

size scaling behaviour of the excess specific heat is governed by h"{x), which 
t^- ■ is proportional to the scaling function /2 discussed in the literature. We 

compare our results with previous work, where the Casimir force has been 
computed by taking the derivative of the excess free energy with respect to 
the thickness of the film. As a preparative study we have also computed the 
scaling functions h!{x) and h{x) for finite L 3 systems with periodic boundary 
conditions in all directions, where L is the linear extension of the system. 
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1 Introduction 



In 1978 Fisher and de Gennes [Tj realized that when thermal fluctuations are re- 
stricted by a container a force acts on the walls of the container. Since this effect 
is rather similar to the Casimir effect, where the restriction of quantum fluctua- 
tions induces a force, it is called "thermodynamic" Casimir effect. Since thermal 
fluctuations only extend to large scales in the neighbourhood of a continuous phase 
transitions it is also called "critical" Casimir effect. Recently this effect has at- 
tracted much attention, since it could be verified for various experimental systems 
and quantitative predictions could be obtained from Monte Carlo simulations of 
spin models [2]. 

In the thermodynamic limit of the three dimensional system, the correlation 
length, which measures the spatial extension of fluctuations, diverges following the 
power law 

£^Q,±i*r (i) 

where t = (T — T c )/T c is the reduced temperature and T c the critical temperature. 
£o,+ and £o,- are the amplitudes of the correlation length in the high and low 
temperature phase, respectively. While £o,+ and £o,- depend on the microscopic 
details of the system, the critical exponent v and the ratio £o,+/£o,- are universal. 
At the critical point also other quantities like the specific heat show a singular 
behaviour: 

C~A±\t\- a + B . (2) 

In the case of the XY universality class that we consider here, the exponent a = 
—0.151(3) [3j of the specific heat is negative. Therefore the analytic background B 
has to be taken into account. Note that the critical exponents of the correlation 
length and the specific heat are related by the hyperscaling relation a = 2 — du, 
where d is the dimension of the system. For reviews on critical phenomena and its 
modern theory, the Renormalization Group, see e.g. [H O EJ [7] . 

The singular behaviour at the critical point originates from the fact that thermal 
fluctuations range over all length scales. Therefore the behaviour in the neighbour- 
hood of the critical point is modified if the system is confined by a container. A 
priori thermodynamic quantities are functions of the reduced temperature and the 
size Lq of the container, assuming a fixed geometry. However the theory of finite 
size scaling predicts that the physics of the system is governed by the ratio L /£ 
as long as L , £ 3> a, where a is the microscopic scale of the system. In particular if 

1 For a review on finite size scaling see [8]. 
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a quantity in the thermodynamic limit behaves as A ~ ao,±|t| w , finite size scaling 
predicts that A(L ,t) ~ LQ^ u g(L /C,), where w is the critical exponent of A and £ 
the correlation length of the bulk system. We can rewrite this equation as 

A(L ,t)^L: /u g(t[L /^n (3) 

by using (CQ), which is the form used in the following. Note that the function g 
depends on the details of the container. For example for a cube it is different from 
a thin film. It also depends on the type of boundary conditions that is imposed by 
the walls of the container on the order parameter of the system. 

The predictions of finite size scaling theory have been tested in experiments 
and theoretical studies for various universality classes and confining geometries; for 
reviews see [8j [9]. Here we shall focus on thin films in the three dimensional XY 
universality class, which is shared by the A-transition of 4 He. Very precise experi- 
mental results for critical exponents and universal amplitude ratios were obtained 
for this phase transition [10] . Also a large number of experiments on thin films 
of 4 He and 3 He- 4 He mixtures were performed to probe finite size scaling [UJ. In 
particular the specific heat C of thin films has been studied. The excess specific 
heat should behave as 

C bulk (t) ~ C(L ,t) ^ L a Q /v ' f 2 {t[Lo/U 1/u ) ■ (4) 

The reason to study the excess specific heat rather than just the specific heat 
C(L Q ,t) is to cancel the analytic background B. Note that the scaling function 
fz(x) of the excess specific heat is, up to a constant factor, the second derivative 
h"{x) of the scaling function h(x) of the excess free energy per area 

f ex = ffiim(L , t) - L f bulk (t) ~ k B TL d+1 h(x) (5) 

where ffii m (L Q , t) is the free energy per area of the thin film, fbuikif) the free energy 
density of the bulk system, d = 3 the dimension of the system and x = t[L /£o] 1/V - 
Note that in the case of thin films we consider, following the literature on the 
thermodynamic Casimir effect, free energies per area. We hope that this does not 
lead to confusion, since in the case of the specific heat, energies per volume are 
considered. 

From a thermodynamic point of view, the Casimir force per unit area is given 

by 

rn df ex , 



where Lq is the thickness of the film. Inserting the finite size scaling ansatz ([5]) for 
the excess free energy into (jBJ) we get 
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This relation is well known and can be found e.g. in [14] . We like to emphasis, that 
it is at the very heart of finite size scaling that the behaviour of the Casimir force, 
which gives the reaction of the film with respect to a change of the thickness and 
the excess specific heat which gives the reaction of the film with respect to a change 
of the temperature are given by the same scaling function h(x) of the excess free 
energy. 

The purpose of the present work is to compute the finite size scaling function 
8(x) by using the energy density of thin films obtained from Monte Carlo simulations 
of a lattice model [T3] and by using experimental data for the specific heat of thin 
films of 4 He near the A-transition [161 E|- 

As a preliminary study, we simulate L 3 systems with periodic boundary condi- 
tions in all directions. In order to eliminate leading corrections to scaling we have 
used the improved two component 4 model on the simple cubic lattice. For a pre- 
cise definition of the model see section [2J Using the results obtained for the energy 
density for a dense grid of temperatures in the neighbourhood of the critical point, 
we investigate the scaling behaviour and compute the finite size scaling functions 
h!{x) and h(x). We find that corrections to scaling are small for the lattices sizes 
F = 8, 16 and 32 that we have simulated. 

In the case of thin films we analyse data for the energy density that were ob- 
tained in [12] from simulations of the improved two component 4 model on the 
simple cubic lattice. In [15] these data where used to compute the specific heat. 
In order to get a vanishing order parameter as it is observed at the boundaries 
of 4 He films, Dirichlet boundary conditions with vanishing field were imposed. In 
singular quantities these lead to corrections oc Lg 1 [TS], which can be expressed 
by an effective thickness £o,e// = L + L s . In [19] we find L s = 1.02(7) for the 
model that we consider here. Note that the boundary conditions also effect the 
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analytic background of the specific heat and the energy density, which also leads 
to corrections oc Lq 1 . However it turns out that these corrections are not given by 
the same Lo, e ff as for the singular quantities. Taking into account these subtleties 
we arrive at accurate result for h'(x), h(x) and 6(x). In particular for 6(x) in the 
range we find a good match with our previous result [20] 0, where we 

computed the Casimir force by taking the derivative of the excess free energy with 
respect to the thickness L of the film. 

Next we compute 8(x) by using experimental results for the excess specific heat 
obtained from experiments on thin films of 4 He [T6t [TT] . Even though this is a quite 
simple exercise, to our knowledge, it has not been done before. For —5 J$ x J$ 4 
we find a reasonable match with our result [20]. However in the low temperature 
phase, for x ;$ — 5 we get results that strongly deviate from [20] and can be ruled 
out by plausibility. This corroborates the observation that in the low temperature 
phase for x ;$ — 5 the excess specific heat does not scale well [TT] . 

This paper is organized as follows: First we define the model and the observables 
that we consider. Next we discuss the finite size scaling behaviour of the free energy 
density. In particular, we discuss corrections to scaling caused by Dirichlet boundary 
conditions. In section H] we compute the scaling functions h'{x) and h(x) for L 3 
systems with periodic boundary conditions. Then in section [5] we compute h'(x), 
h(x) and 8(x) using the data for the energy density of thin films with Dirichlet 
boundary conditions obtained in [13] . The result for 8(x) is compared with the 
one that we [20] obtained directly from the thermodynamic Casimir force. Next in 
section [6] we compute 9{x) starting from data for the excess specific heat of films of 
4 He in the neighbourhood of the A-transition. Finally we summarize and conclude. 

2 The model and the observables 

We study the two component 4 model on the simple cubic lattice. We label the sites 
of the lattice by x = (x ,Xi,x 2 ). The components of x might assume the values 
Xi G {1, 2, . . . , Li}. In this work we have performed simulations of lattices with 
L = Li = L 2 and periodic boundary conditions in all three directions. Furthermore 
we analyse data obtained in [15] for thin films. In this case lattices of the size 
Li = L 2 = L and L <C L are studied. In 1 and 2-direction periodic boundary 
conditions and free boundary conditions in 0-direction are employed. This means 

2 In [20] we have compared our result for 0(x) with previous ones obtained from simulations of 
the XY model [211 E2] and experiments on thin films of 4 He [TJl [13] ; overall we find a reasonable 
agreement. 
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that the sites with xq = 1 and xq = Lq have only five nearest neighbours. This type 
of boundary conditions could be interpreted as Dirichlet boundary conditions with 
as value of the field at xq = and xq — Lq + 1. Note that viewed this way, the 
thickness of the film is Lq + 1 rather than Lq. This provides a natural explanation 
of the result L s = 1.02(7) obtained in [19] . The Hamiltonian of the two component 
4 model, for a vanishing external field, is given by 



where the field variable <p x is a vector with two real components. < x, y > denotes 
a pair of nearest neighbour sites on the lattice. The partition function is given by 



Note that following the conventions of our previous work, e.g. [23] , we have absorbed 
the inverse temperature (3 into the Hamiltonian. [f| In the limit A — ► oo the field 
variables are fixed to unit length; i.e. the XY model is recovered. For A = we get 
the exactly solvable Gaussian model. For < A < oo the model undergoes a second 
order phase transition that belongs to the XY universality class. Numerically, 
using Monte Carlo simulations and high-temperature series expansions, it has been 
shown that there is a value A* > 0, where leading corrections to scaling vanish. 
Numerical estimates of A* given in the literature are A* = 2.10(6) [21], A* = 2.07(5) 
[23J and most recently A* = 2.15(5) [3]. The inverse of the critical temperature 
f3 c has been determined accurately for several values of A using finite size scaling 
(FSS) [3]. We shall perform our simulations at A = 2.1, since for this value of A 
comprehensive Monte Carlo studies of the three-dimensional system in the low and 
the high temperature phase have been performed [191 El ESI EE]. At A = 2.1 one 
gets p c = 0.5091503(6) [3J. Since A = 2.1 is not exactly equal to A*, there are still 
corrections oc L~~ w , although with a small amplitude. In fact, following [3], it should 
be by at least a factor 20 smaller than for the standard XY model. 

In [TJ] we find for A = 2.1 by fitting the data for the second moment cor- 
relation length in the high temperature phase ^2nd — 0.26362(8)t~ ' 6717 , where 
t = 0.5091503 — (3. We shall use this definition of the reduced temperature also 
in the following discussion of our numerical results; Hence £o = 0.26362(8). Note 
that in the high temperature phase there is little difference between £ 2n d and the 

3 Therefore, following [6] we actually should call it reduced Hamiltonian. 
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exponential correlation length £ exp which is defined by the asymptotic decay of the 
two-point correlation function. Following [23] lim^o |^ — 1-000204(3) for the 
thermodynamic limit of the three-dimensional system. Hence at the level of preci- 
sion reached here it does not matter whether ^Q >exp or £o,2nd is used in the scaling 
variable x = t[L /£o] 1/iy - 



2.1 The internal energy and the free energy 

The reduced free energy density is defined as 

/ = -— S-p- logZ. (11) 

I.e. compared with the free energy density /, a factor fc#T is skipped. 

Note that in eq. (Q (3 does not multiply the second term. Therefore, strictly 
speaking, (3 is not the inverse of k^T . In order to study universal quantities it is 
not crucial how the transition line in the f3-\ plane is crossed, as long as this path is 
not tangent to the transition line. Therefore, following computational convenience, 
we vary (3 at fixed A. Correspondingly we define the (internal) energy density as 
the derivative of the reduced free energy density with respect to (3. Furthermore, 
to be consistent with our previous work [15], we multiply by —1: 

E = 1 *°« Z . (12) 
L LiL 2 dp v ; 

It follows 



(13) 

which can be easily determined in Monte Carlo simulations. From eqs. ( fl~TlfT2l) it 
follows that the free energy density can be computed as 

f(P) = fW- f ■ (14) 

J 00 



3 The finite size scaling behaviour of the free en- 
ergy 

Let us briefly discuss the scaling behaviour of the reduced excess free energy per 
area. Since we study an improved model we ignore corrections oc Lq U in the follow- 
ing. We take into account leading corrections due to the boundary conditions by 



6 



replacing the thickness L of the film by £o,e// = L + L s at the appropriate places. 
We split the free energies in singular (s) and non-singular (ns) parts: 

fex{t,L ) = ffU m (t, Lq) — Loffaikit) 

— ffilm,s(t, Lq) + Lo,eff,nsfns(t) ~ ^oiWfc.sW — Lof ns (t) 

= L~ 2 eff h(x) + L s (t)f bu i ktS (t) + L sns (t)f ns (t) (15) 

where 

h(x) = L Qe ff[ffii m)S (t,Lo) — Lo, e fffbuik,s(t)\ (16) 

is a universal finite size scaling function and x = t[Lo, e ///£o] 1/ ' ! ''- Following RG 
theory the non-singular part is not affected by finite size effects. However it is not 
clear a priori how Dirichlet boundary conditions affect the non-singular part of the 
free energy. Therefore we allow for L sns = L 0fi jf^ ns — L 7^ and L sns 7^ L s . Taking 
the derivative with respect to L we get the thermodynamic Casimir force per area 

m 

/3F casimir = - ^ — - = 2L ^ ff h(x) — L Q ^ ff -xh'(x) = L Q 3 eff 6(x) (17) 
where 

6{x) = 2h{x) - -xh'{x) . (18) 

Note that the boundary terms L s fb u ik, s and L sns f ns do not contribute to the Casimir 
force. 

4 Warmup exercise: finite cubic system with pe- 
riodic boundary conditions 

First we have studied a finite lattice with L = L = L\ = L 2 with periodic bound- 
ary conditions in all directions. This way we avoid corrections caused by the free 
boundary conditions and possible difficulties related with the Kosterlitz-Thouless 
transition of thin films. In our simulations we determine the energy density E for 
the lattice sizes L — 8, 16 and 32 for a large number of /3-values in the neighbour- 
hood of the critical temperature. In particular, we have simulated at 159 /9-values 
in the interval 0.35 < (3 < 0.58, 205 /3-values in the interval 0.35 < (3 < 0.58 and 
161 /3-values in the interval 0.45 < (3 < 0.535 in the case of L — 8, 16 and 32, 
respectively. For most of these simulations we performed 10 6 measurements. For 
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each of these measurements we performed one Metropolis sweep, two overrelaxation 
sweeps and a number of single cluster [27] updates. The number of single cluster 
updates is chosen such that the number of updates times the average size of a cluster 
roughly equals L 3 . As random number generator we have used the SIMD-oriented 
Fast Mersenne Twister algorithm [28]. In total we have used about 3 weeks of 
CPU-time on a single core of a Quad-Core Opteron(tm) 2378 CPU (2.4 GHz). In 
order to compute the excess energy density 

E ex (L,t) = E(L,t) — E bu i k (t) (19) 

we have used the results for E bu ik(t) obtained in section 4.1 of [15] . 
The reduced excess free energy density behaves as 

f(t,L)-f bulk (t)~L- 3 h(x) (20) 

where x = t[L/^Q\ l / u . It follows for the excess energy density 

E(t, L) - E bulk (t) ~ [L/^L- 3 h'(x) . (21) 

Since we study an improved model and periodic boundary conditions in all di- 
rections, we expect that analytic corrections are leading. In figure [T] we plot 
[L/£ r 1/ly £ 3 [£(£,£) - E bulk (t)} as a function of t[L/^ ]-^ u . The curves for L = 8, 
16 and 32 fall nicely on top of each other, showing that corrections to scaling are 
numerically small. The excess energy is positive for all temperatures. At x ~ —0.3 
the function assumes a maximum. The value at the maximum is ~ 0.295. In the 
high temperature phase, for increasing x the function h'(x) rapidly approaches zero. 
In contrast, in the low temperature phase it is only slowly approaching zero with 
decreasing x. This behaviour might be explained by the presence of a Goldstone 
mode in the low temperature phase. 

Next we have computed the excess free energy using eq. (1141) . To this end we 
integrated our data for the excess energy by using the trapezoidal rule. We have 
started the integration at (3q = 0.35, 0.45 and 0.49 for L — 8, 16 and 32, respectively. 
At these values of (3q the deviation of the excess energy from zero is of similar size 
as the statistical error. In figure [2] we plot L 3 [f(t,L) — f bu ik(t)} as a function of 
t[L/£o] -1 /^ '. As one should expect, the curves for L = 8, 16 and 32 fall nicely on 
top of each other. Since the excess energy is positive for all values of x, the excess 
free energy is monotonically increasing with increasing x. At the critical point the 
scaling functions assumes the value h(0) = —1.162(4). 

In Fig. 6 of [29] results for an isotropic L 3 system with periodic boundary 
conditions for the Ising universality class and the limit n —>■ oo of 0(n) symmetric 
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Figure 1: We consider an L 3 system with periodic boundary conditions in all 
directions. We plot [E(t,L) - E hM (t)][L/£ Q ]- 1 / v L 3 as a function of t(L/£ ) 1/y for 
L — 8, 16 and 32, where we use v = 0.6717 and £ = 0.26362. For a discussion see 
the text. 
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Figure 2: We consider an L 3 system with periodic boundary conditions in all 
directions. We plot L 3 [f(t,L) — fbuik{t)} as a function of t(L /^ ) 1 ^ 1 ' for L = 8, 16 
and 32, where we use v = 0.6717 and £ = 0.26362. For a discussion see the text. 
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models are given. The Ising result is obtained from a perturbative approach in three 
dimensions fixed, while the n — > oo result is exact. In the Ising case, h(x) shows 
a minimum in the low temperature phase, while for n — > oo it is monotonically 
decreasing as the temperature decreases. Hence qualitatively the behaviour for 
n = 2, studied here, is the same as for n — > oo. 

5 Film geometry with free boundary conditions 

Here we study thin films with free boundary conditions in the short direction. This 
geometry is relevant for the comparison with experimental results obtained for thin 
films of 4 He. Most of the Monte Carlo data are taken from our previous work [T5J [20] , 
where we have simulated films of the thicknesses L = 8, 16 and 32. Therefore we 
refrain from giving the details of the simulations and refer the reader to [H2, [20] . 

Analogous to the previous section we compute the scaling function h'(x) of the 
excess energy density and h(x) of the excess free energy density. Using these we 
obtain the scaling function 8(x) = 2h(x) — f;h'(x) of the thermodynamic Casimir 
force. 

In [15] we have taken great care to get the deviations of the energy density 
from its effectively two dimensional thermodynamic limit under control. In order 
to achieve this, quite large ratios L\/Lq are needed in the neighbourhood of the 
peak of the specific heat. In the case of L = 8 we have simulated lattices of a size 
up to L\ = L 2 = 2048, and for L Q = 16 up to Li = L 2 = 1800. For L = 32 we have 
skipped the interval 0.5136 < (3 < 0.516, since we could not simulate sufficiently 
large values of L\ = L 2 . 

In figure [3], similar to the previous section, we have plotted E ex Ll[L / '^o] -1 ^ 
versus t[L /^o] 1//!/ , where now 

E ex (L ,t) = L [E(L ,t) — Ebuik(t)] (22) 

is the excess energy per area. In contrast to the previous section we find that there 
is a huge discrepancy between the three curves. The dominant effect seems that the 
curves are shifted by a constant with respect to each other. Note that replacing Lq 
by Lo } eff = L + L s with L s = 1.02(7) [19] does change this situation only little. In 
section [3] we have argued that the analytic background of the energy density might 
suffer from a boundary correction that is not given by the effective thickness -^o,e// — 
L +L s that describes the leading boundary corrections of singular quantities. Below 
we shall study this question in detail at the critical point of the bulk system, where 
we have data for thicknesses up to Lq = 64 available. 
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Figure 3: We plot E ex Ll[L /£ ] x l v as a function of ijXo/^o] 1 ^ f° r thin films, where 
we have used £ = 0.26362 and v = 0.6717. For a discussion see the text. 
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Table 1: Energy density of thin Rims of the thickness Lo at the inverse critical 
temperature (3 C = 0.5091503(6) of the three dimensional systems. In all cases 
L\ — L2 = 6Lq. 





E 


8 


0.799566(31) 


12 


0.832786(19) 


16 


0.850727(13) 


24 


0.8698028(85) 


32 


0.8798552(57) 


48 


0.8903321(37) 


64 


0.8957662(29) 



5.1 Finite Size Scaling at the critical point of the bulk sys- 
tem 

In order to get a better understanding of the corrections we have studied in detail 
the behaviour at the critical point of the three dimensional bulk system. In the 
context of [15] we have simulated lattices of the thickness Lo = 8, 12, 16, 24, 32, 
48 and 64 and L\ = L2 = QLq. In [15J we have checked that this choice of L%, L2 is 
sufficient to approximate well the effectively two dimensional thermodynamic limit 
of the film at the critical point of the three dimensional system. Our results for the 
energy density are summarized in table [TJ 

In the case of periodic boundary conditions in all directions, the energy density 
at the critical point behaves as 

E(L) = E ns + cL~ d+1 ^ (23) 

where d = 3 is the dimension of the system. Using lattices of the size Lo = L\ = L2 
with L up to 128 we find [25] 

E ns = 0.913213(5) + 20 x (J3 C - 0.5091503) + 5 x 10~ 7 x (1/a + 1/0.0151) . (24) 

In order to take into account corrections due to the free boundary conditions of 
the thin film we use the ansatz 

L E(L ) = (L + L sns )E ns + c f L- d +} +1/ " (25) 
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to fit the data given in tabled! As input we have used v = 0.6717(1) [3], E ns given 
in eq. (124")) and L s = 1.02(7) [19]. The parameters of the fit are L sns and cj. Fitting 
all data with Lq > 8 we get an acceptable xV d -°- f - We find L ™ = -1.3529(3) 
when fixing L s = 1.02 and L ns = —1.3523(3) fixing L s = 0.95. Hence L sns is clearly 
different from L s and it shows little dependence on the value taken for L s . We have 
checked that the error of L sns due to the uncertainties of v and E ns can be ignored. 

5.2 Taking into account boundary corrections 

The boundary correction L ns E ns should be an analytic function of the reduced 
temperature. In a first attempt we shall approximate it by its value at the critical 
point of the three dimensional system found above. Hence in figure H] we plot 

^e X Ll eff [L 0teff /^ ]- 1 / v where 

E ex = L Q E(L Q ,t) — L ^ e ffE bu i k (t) + {L s — L sns )E ns (26) 

as a function of x = tfLo^/y/^o] 1 ^- Now we see a quite good matching of the three 
curves. Only for small x discrepancies are visible. 

Next we have computed the finite size scaling function 9 of the Casimir force 
following eq. ( fi~8l) . In the case of Lq = 8 and Lq = 16 we have used the function 
h'(x) as given in figure HI In the case of Lq = 32 we have taken the missing part 
in the range —5.9 > x > —9.1 from the results for Lq = 16. To this end we 
have matched the values of the function at x — —5.9 and x = —9.1 resulting in 
h' 32 (x) = h' 16 (x) + 0.011 - 0.002(a; + 5.9) for -5.9 > x > -9.1. We have computed 
the function h(x) by numerically integrating h'(x) using the trapezoidal rule. For 
sufficiently large x the Casimir force vanishes and therefore h(x) = ^h'(x). Hence 
for large x: 

h\x) = cx 2u ~ l . (27) 

In [20] we found that the thermodynamic Casimir force is of similar size or smaller 
than the numerical errors that we achieve for x ^ 4. We have checked that in this 
range the scaling function h'(x) of the excess energy indeed follows eq. (I2TI) . 

Hence we have started the numerical integration in the high temperature phase 
at x ~ 4 with the starting value h(x ) = ^h'(x ). In order to check the reliability 
of our result, we have redone the integration using a set of data points, where we 
have skipped every second value of (3. We found an agreement within the statistical 
errors. Our results for 6 are plotted in figure El In the range — 7 < x < 5 the 
curves obtained from the different values of L match quite well. There is also a 
good match with 9 obtained in [20]. In particular the value and the position of 
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Figure 4: We plot E ex Ll eff [L 0te ff /£ ] x l v as a function of t[L 0je ff / '^o] 1 ^ for L = 8, 
16 and 32. For a discussion see the text. 
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the minimum of 6 are completely consistent. However for x < — 7 the difference 
between the curves becomes clearly visible and increases with decreasing x. For 
small x, even for Lq = 32 there is a huge discrepancy with the result of [20J. 

Since these discrepancies appear for large values of |x| it is likely that they are 
mainly caused by analytic corrections. To check this explicitly, we allowed for two 
different types of corrections: 

x = t(l-ct)[L , eff /^] 1/u (28) 

and for a temperature dependence of the boundary correction of the analytic part 
of the energy 

E ex = L E(L , t) — Lo te ffEb u ik(t) + (L s — L sns )E ns — c^t . (29) 

We find that the curves for h'(x) obtained from Lq = 16 and Lq = 32 can be 
nicely matched by adjusting the two parameters c and c&. Matching in the interval 

— 18 < x < 3 we find c ~ —1.1 and c& m —3.03 and for the interval —25 < x < 5 
c ~ —0.75 and q, ~ —2.97. Using the corresponding results for h'(x) we have 
computed the finite size scaling function 9(x) that is plotted in figure [6j Now 
we see that the rage of the matching with our previous result |20j is extended to 

— 15 <; x <J 4 in the case of the matching range —25 < x < 4. For still smaller 
values of x discrepancies rapidly increase. Likely higher order analytic corrections 
are the main reason for this behaviour. However also other types of corrections like 
t uw with uj' 1.8 [30] should be taken into account. Therefore we abstain from 
fitting with t 2 corrections. 

6 The specific heat of thin films of 4 He and the 
thermodynamic Casimir force 

In a number of experiments the excess specific heat of thin films of 4 He and 3 He- 4 He 
mixtures has been measured in the neighbourhood of the A-transition [llj . In these 
works the scaling function /*2 which is defined by 

CU(f) - C{L , t) ~ L% /u f 2 (tLl /u ) (30) 

is extracted from experimental data for the specific heat of the three-dimensional 
bulk system Cb u ik{t) and of thin films C(Lq, t), where Lq is the thickness of the film. 
Since the specific heat is the derivative of the energy density with respect to the 
temperature, fzix) is, up to a constant factor, equal to h"{x). 
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Figure 5: We plot finite size scaling function 9 of the thermodynamic Casimir force. 
We have computed 9 from the finite size scaling function h' of the excess energy 
per area of the film. For comparison we give the result of our previous work [20J 
where we have computed 9 directly from the thermodynamic Casimir force. For a 
discussion see the text. 
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Figure 6: Similar to the previous figure. Here we have taken into account analytic 
correction computing the scaling function h! . For a discussion see the text. 
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To compute this factor let us start from the excess reduced free energy density: 

f(t,L )-f bulk (t)~L 3 h(x) (31) 

where t = T/T\ — 1 and x = tfLo/^o] 1 ^- Note that here, as long as the free energy 
density and Lq 3 are measured in the same units, h(x) is uniquely defined; there is 
no ambiguous factor. 

The excess energy density is given by the derivative with respect to (3 = l/(k B T). 
Hence 

E(t, L ) - E bulk {t) = -L 3 [L /^-^—(3- 2 ~ -L 3 [L /^] 1/u k B T x (32) 

where we have approximated (3~ 2 ~ &b^a in the neighbourhood of the A-transition. 
The specific heat as defined in the experiments is given by the derivative of the 
energy density with respect to the temperature. Hence 

C bulk {t) - C(t, Lo) ~ k B L^[L Q /^h"{x) . (33) 

The results for the specific heat of the experiment are given in JmoT l K- 1 . 
These we convert into k B A~ 3 to get the same units on both sides of equation (|33|) . 
Note that the thickness of the films in [T6l[T7] is quoted in A -3 . To this end we need 
the density p\ = 146.1087 kg/m 3 [31J of 4 He at the A-transition, the molar weight 
4.0026 . . . g mol™ 1 of 4 He and the Boltzmann constant k B = 1.38065 . . . x KT 23 JK _1 . 
This amounts to the factor 0.00264 . . . J~ l mo\Kk B A~ 3 . In [32J the data are given 
as a function of the reduced temperature t — 1 — T/T\, where T A = 2.17 . . . K . In 
order to plot them function of x = t{L /£ ] 1/u we have used £ = 1.422(5) A 
which we [15] have computed from the amplitude of the bulk specific heat of 4 He 
at vapour pressure [33] and the universal amplitude ratio R a [26|. 

In figure [7] we have plotted 

[C bulk (t) - C(L ,t)}k B l L 3 [L /^ Q ]- 2 ^ (34) 

as a function of x = t[L / '^o] 1 ^ . To this end we have used the data given in [32J for 
the thicknesses 483, 1074, 2113, 5039, 6918 and 9869A. 

Note that h"(0) ~ 3.05 as can be obtained from the results of section 4.3 of [15J. 
For x <J —5 the curves obtained from different thicknesses of the film fall reasonably 
well on top of each other. It has been noticed [H] that for x J$ —5, in particular in 

4 It would be interesting to repeat the analysis for other data sets as e.g. those of [33]. 
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Figure 7: We plot - C(L , t)][L /^ ]- a/u as a function of x = t[L /^] 1/u - 

The data for thin films of 4 He of thicknesses 483, 1074, 2113, 5039, 6918 and 9869A 
obtained in [T6l [T7] are taken from [32] . Note that at the critical point t — the 
finite size scaling function assumes the value ~ 3.05. For a discussion see the text. 
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the neighbourhood of the minimum, the results obtained for different thicknesses 
differ by quite large factors. In [15] we have computed the specific heat and the 
scaling function /2 starting from the data for the energy density discussed above. 
We find that for x J$ —5 our final result is clearly smaller than the experimental 
ones [161CE7]. 

Starting from the results for the finite size scaling function h"(x) obtained from 
different thicknesses we have computed the scaling function h'(x). To this end, 
we have applied the trapezoidal rule. Similar to the previous section, we have 
started the integration at xq ~ 4 in the high temperature phase. As starting value 
we have taken h'(xo) = -^—^h"{xo). Again we have integrated h'(x) using the 
trapezoidal rule to get h(x). Here we have taken the same value for xq as above and 
h(xo) = 7^h'(xo). Using these results for h\x) and h(x) we have computed 9(x) 
which we have plotted in figure [HJ 

In order to check the effect of errors due to the finite step-size of the integration, 
we have repeated the integration, skipping every second value of x. In order to 
check the effect of the singularity of h'(x) we have fitted the data for the specific 
heat in the neighbourhood of the transition with the ansatz 

h"(x) =3.Q5 + c±\x\- a (35) 

Then we have integrated the ansatz in the neighbourhood of the transition and 
compared it with the corresponding result from the trapezoidal rule. We find that 
the numerical results only change little and the conclusions drawn below are not 
effected. 

Let us now discuss the results that we have obtained: For 483A the curve 
is monotonically decreasing with decreasing x; in particular no minimum of the 
function can be observed. For 1074A a shallow minimum occurs at x ~ —5.3; for 
x ^ 6.1 the function is decreasing again with decreasing x. For 2113A we see a clear 
minimum at x ~ —5.5; However the value of the minimum is clearly smaller than 
the one of [20J. In the case of 5039 A and 69 18 A we find a quite good match with 
our result [20] down to x ~ —7. For 6918A the minimum is located at x ~ —4.8 
and the value of the minimum is 9 « —1.3. For 6918A the minimum occurs at 
x ~ 4.65 and the value of the minimum is 9 « —1.36. For 5039A no data for 
x < —7 are available. For 6918A the curve is decreasing again for x J$ —6.7 with 
decreasing x. Up to here the expectation that with increasing thickness of the 
film the result converges toward the universal scaling function is fulfilled. However 
for the largest thicknesses studied, 9869A even the worst mismatch is found. The 
curve is monotonically decreasing with decreasing x and even for x ^ 5 there is 
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Figure 8: We plot results for the scaling function 9 of the thermal Casimir force. 
These results were obtained by numerically integrating experimental data [TH [T7J 
[32] for the excess specific heat of thin films of 4 He near the A-transition. For 
comparison we give the result obtained in [20J. For a discussion see the text. 
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quite large mismatch with our result for 9 [20]. Playing around with the data, 
one finds that smaller values of h"(x) for needed to avoid that 9 is 

decreasing with decreasing x for x J$ —5. One should note that eq. (El) only holds 
for the scaling limit. Therefore the observations made here are not an indication 
that the experimental data are effected by an error. They can also be explained by 
corrections to the scaling behaviour. This is at least true for the smaller thicknesses 
like 483A. The result for 9869A is more puzzling. One should note that the analysis 
presented in this section does not depend on the type of boundary conditions that 
is realized in the experiment. But we think that the behaviour of 9{x) for x J$ —5 
obtained here can not be explained by different boundary conditions from those 
used in the study of the lattice model [151 120] . 

7 Summary and Conclusions 

We have studied the relation of the excess specific heat, the excess energy and the 
thermodynamic Casimir force in thin films in the three-dimensional XY universality 
class. To this end we have exploited the relation (|Sj) 

9{x) = 2h(x) - -ti{x) (36) 

among the finite size scaling functions 9(x) of the thermodynamic Casimir force, 
h(x) of the excess free energy and h\x) of the excess energy. We have analysed 
data obtained for the energy density of the improved two-component 4 model on 
the simple cubic lattice and experimental results for the specific heat of thin films 
of 4 He near the A-transition [161 E] • 

As a first exercise we have computed the functions h'(x) and h(x) for a finite 
lattice of the size L 3 with periodic boundary conditions in all directions. This way 
we avoid potential problems related with the Kosterlitz-Thouless phase transition 
of the thin film and corrections caused by Dirichlet boundary conditions. Indeed 
we find a good collapse of the data already for rather small lattice sizes L — 8, 16 
and 32. 

Next we have repeated the same exercise for films with free boundary conditions 
using the data for the excess energy density obtained in [15] . Here we find a huge 
mismatch between the thicknesses L = 8, 16 and 32. Replacing L by Lo,e// = 
L + L s with L s = 1.02(7) [19] does not remove this discrepancy. We argue that 
the non-singular part of the energy density suffers from boundary corrections that 
are not described by the same L s which accounts for the corrections in singular 
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quantities. Indeed, the analysis of the energy density of films up to the thickness 
L = 64 at the critical temperature of the three dimensional system results in a 
L sns for the analytic background of the energy density that is clearly different from 
L s . Taking into account this result we find a reasonable collapse of the functions 
obtained from Lq = 8, 16 and 32 in a large range of the scaling variable x. 

Computing the Casimir force from this result for h'(x) we find a good collapse 
for —7 J$ x < 4. In this range of x we also find a good agreement with the result for 
9 that we have obtained in [20] . Note that in particular the minimum of the scaling 
function 9 is within this range. We confirm the position x m i n and the value 9 m i n 
that we have obtained in [20J. Next we took into account analytic corrects. The 
coefficients of these corrections were computed by matching the results obtained 
from L = 16 and L = 32. This way we could extend the range of agreement with 
our previous result |2pJ down to x ~ —15. 

We have presented a viable alternative to compute the scaling function 9 of the 
Casimir force. The nice agreement with our previous result [2U] gives us further 
confidence in the correctness of the results. 

Finally we have computed the scaling function 9 using experimental results for 
the excess specific heat [IB], [TT] . Here we find a reasonable match with the theoretical 
results in the range x ^ —5. In particular we find evidence that the minimum of the 
scaling function 9 is located at x m i n ~ 5 which is consistent with our prediction [20] 
but slightly larger than x m i n = —5.45(12) p2] and x m i n = —5.7(5) [13J, where the 
thermodynamic Casimir force has been determined for films of 4 He with thicknesses 
< 600A. 

In the range x J$ —5 the curves obtained from different thicknesses show quite 
different behaviour. Furthermore in this range 9(x) is decreasing with decreasing 
x. This behaviour corresponds to too large values of h"(x) in the range x J$ —5. 
We think that understanding this problem requires a detailed knowledge of the 
experiments and is therefore beyond the scope of the present work. 
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